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5£ ■ Abstract 

We used an Artificial Neural Network (ANN) to derive the orbital parameters of spectro- 
scopic binary stars. Using measured radial velocity data of seven double-lined spectroscopic 
binary systems V373 Cas, V2388 Oph, V401 Cyg, GM Dra, V523 Cas, AB And, and HD 
141929, we found corresponding orbital and spectroscopic elements. Our numerical results 
are in good agreement with those obtained by others using more traditional methods. 
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^ ■ 1 Introduction 



Determining the orbital elements of binary stars helps us to obtain fundamental information, 
^ ■ such as the masses and radii of individual stars, that has an important role in understanding 

the present state and evolution of many interesting stellar objects. Analyzing of both light and 
radial velocity (hereafter RV) curves, derived from photometric and spectroscopic observations, 
respectively, yields a complete set of basic absolute parameters. 

There are different methods to determine the orbit of a spectroscopic binary from its RV 
curve. Lehmann-Filhes (1894) introduced a geometrical method to determine the orbital ele- 
ments from the geometrical properties of the RV curve, especially its maxima and minima. The 
method of Lehmann-Filhes has been found to be very useful, and little, if any, longer than other 
methods, providing a planimeter is used. This method also gives, for improving the final solu- 
tion of the orbit, the differential corrections to the preliminary elements using the form of the 
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equations of condition obtained by the method of least squares (see also Petrie 1960). Sterne 
(1941) described two forms of least square solutions. The first form is suitable for all orbits 
except those with very small eccentricities. The second form is particularly suitable for orbits 
having very small eccentricities (see also Petrie 1960). Karami & Teimoorinia (2007) introduced 
a new non-linear least squares velocity curve analysis technique for spectroscopic binary stars. 
Their method was applicable to orbits of all eccentricities and inclination angles and the time 
consumed was considerably less than the method of Lehmann-Filhes. They showed the validity 
of their new method to a wide range of different types of binary. See Karami & Mohebi (2007a, b) 
and Karami et al. (2008). 

In the present paper, we use an Artificial Neural Network (ANN) to find the optimum match 
to the four parameters of the RV curves of the seven double-lined spectroscopic binary systems: 
V373 Cas, V2388 Oph, V401 Cyg, GM Dra, V523 Cas, AB And, and HD 141929. Our aim is 
to show the validity of our new method for a wide range of different binary types. 

The spectral type of the primary and secondary component of V373 Cas is B0. 511 and B4III, 
respectively. The mean effective temperature is 22000 K and 18000 K for the primary and 
secondary components. The angle of inclination is ~ 60° with a period of 13.4 days (Hill & 
Fisher 1987). V2388 Oph is very close visual binary. The spectral type is F3 V with a relatively 
long period of 0.802 days. The system appears to be one of the most luminous among currently 
known contact binaries. Orbital inclination angle is 90° (Rucinski et al. 2002). V401 Cyg 
appears to be a rather typical contact system. The orbital period is 0.582714 days (Rucinski 
et al. 2002). GM Dra is a contact system of the W-type, with spectral type F5V and period 
0.338741 days (Rucinski et al. 2002). V523 Cas is one of the faintest known contact binaries. 
The spectral type is K4V and the period is 0.233693 days (Rucinski et al. 2003). AB And is 
a contact binary with spectral type G8V and period 0.3318919 days (Pych et al. 2004). HD 
141929 is a double-lined spectroscopic binary with a period of 49.699 days. The orbit is eccentric 
(e = 0.393). The effective temperature for both components is estimated to 9500 ± 250 K. Both 
components have the same spectral type A0/1V, however, the secondary is rotating slower than 
the primary and the inclination of the orbit is about 11° (Carrier 2002). 

Following Smart (1990), the radial velocity of a star in a binary system is defined as follows 

RV = V cm + K[cos(6 + uj) + e cos u] , (1) 

where V cm is the radial velocity of the center of mass of system with respect to the sun. Where 
6, uj and e are the angular polar coordinate (true anomaly), the longitude of periastron and the 
eccentricity, respectively. Also 

2ir asini . . 

where P is the period of motion, a is the semimajor axis of the orbit and inclination i is the 
angle between the line of sight and the normal of the orbital plane. 

From Smart (1990) for a small eccentricity, e <C 1, the true anomaly, 9, can be expressed in 
terms of the photometric phase, 4>, as follows 

e 3 5 
6 = 2ir<f> + (2e - — ) sin(2vr0) + -e 2 sin(4vr0) 

+ l^ e 3 sin ( 67r 0) + o( e 4 ). (3) 

Here we apply the ANN method to estimate the four orbital parameters, V cm , K, e and uj 
of the RV curve in Eq. (1). ANNs have become a popular tool in almost every field of science. 
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In recent years, ANNs have been widely used in astronomy for applications such as star/galaxy 
discrimination, morphological classification of galaxies, and spectral classification of stars (see 
Bazarghan et al. 2008 and references therein). Following Bazarghan et al. (2008), we employ 
Probabilistic Neural Networks (PNNs). An example of a PNN is shown in Fig. 1. This network 
has been investigated in ample details by Bazarghan et al. (2008). 

In this work, for the identification of the observational RV curves, the input vector in Fig. 1, 
X = (X\,X2, X n ), is the fitted RV curve of a star with 36 data points (n = 36). This number 
is enough to cover the whole of RV curve. The network is first trained to classify RV curves 
corresponding to all the possible combinations of V cm , K, e and lo. For this we synthetically 
generate RV curves given by Eq. (1). We generate one RV curve for each combination of the 
parameters: 

• -100 < V cm < 100 in steps of 1; 

• 1 < K < 300 in steps of 1; 

• < e < 1 in steps of 0.001; 

• < uj < 360° in steps of 5. 

Note that from Petrie (1960), one can guess V cm , K and e from a RV curve. This enable us 
to limit the range of parameters around their initial guesses. This gives a set of pattern groups, 
one group for each combination of V cm , K, e and to. Each pattern group, k, is characterized 
by N/, Gaussian functions (see Bazarghan et al. 2008). For each system, we first fit a curve on 
the observational RV data. Then using the fitted curve, the RV is computed in 36 photometric 
phases. When a observational RV curve of an unknown classification is fed to the network, the 
summation layer of the network computes the probability functions of each class. Finally at 
the output layer we have C, the value with the highest probability (see again Fig. 1). 

2 Numerical Results 

Here PNN is used as a tool to derive the orbital parameters of the seven different double-lined 
spectroscopic systems V373 Cas, V2388 Oph, V401 Cyg, GM Dra, V523 Cas, AB And, and HD 
141929. Using measured RV data of the two components of these systems obtained by Hill & 
Fisher (1987) for V373 Cas, Rucinski et al. (2002) for V2388 Oph, V401 Cyg and GM Dra, 
Rucinski et al. (2003) for V523 Cas, Pych et al. (2004) for AB And, and Carrier (2002) for HD 
141929, the fitted velocity curves are plotted in terms of the photometric phase in Figs. 2 to 8. 

The orbital parameters obtaining from the ANN for V373 Cas, V2388 Oph, V401 Cyg, GM 
Dra, V523 Cas, AB And, and HD 141929 are tabulated in Tables 1, 3, 5, 7, 9, 11 and 13, 
respectively. Tables show that the results are in good accordance with the those obtained by 
aforementioned authors. Tables 9 and 11 show that the results of eccentricities for V523 Cas 
and AB And are not significantly different from zero. A Monte Carlo analysis clears that they 
should really be assigned circular orbits. 

The combined spectroscopic elements including m p sin 3 i , m s sin 3 i, (a p + a s ) sin i and m p / m s 
are calculated by substituting the estimated parameters K, e and lo into Eqs. (3), (15) and 
(16) in Karami & Teimoorinia (2007). The results obtained for the seven systems are tabulated 
in Tables 2, 4, 6, 8, 10, 12 and 14 show that our results are in good agreement with the those 
obtained by aforementioned authors. 

Note that the errors in the tables of orbital parameters and combined spectroscopic elements 
in Karami & Mohebi (2009) seem to be unrealistic errors. But they are indeed the standard 
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errors which are obtained from the nonlinear least squares of Eq. (14) in Karami & Teimoorinia 
(2007). The errors of observational RV data are not included in them. Following Lucy & 
Sweeney (1971) the meaningful errors <jj; i = 1,2,3,4 corresponding to the orbital elements 
V cm , K,e,to, respectively, can be obtained from the inverse of the error diagonal matrix A as 
of = 0-2^-1 with An = N ^ A<22 = A u /2, A 33 = NK 2 /2 and A AA = (7re/180) 2 A 33 . Where 

cr is the standard error of an observation of average weight. Using Eqs. (3), (15) and (16) in 
Karami & Teimoorinia (2007), the errors of combined spectroscopic elements are obtained from 
the errors of orbital parameters. 

3 Conclusions 

An Artificial Neural Network to derive the orbital elements of spectroscopic binary stars is 
applied. This method is applicable to orbits of all eccentricities and inclination angles of different 
types of binaries. In this method the time consumed is considerably less than the method of 
Lehmann-Filhes and even less than the non-linear regression method introduced by Karami & 
Teimoorinia (2007). It is possible to make adjustments in the elements before the final result 
is obtained. There are some cases, for which the geometrical methods are inapplicable, and in 
these cases the present one may be found useful. One such case would occur when observations 
are incomplete because certain phases could not have been observed. Another case in which 
this method is useful is that of a star attended by two dark companions with commensurable 
periods. In this case the resultant velocity curve may have several unequal maxima and the 
geometrical methods fail altogether. 

Using the measured RV data of V373 Cas, V2388 Oph, V401 Cyg, GM Dra, V523 Cas, AB 
And, and HD 141929 given by aforementioned authors, we find the orbital elements of these 
systems by the PNN. Our numerical results shows that the results obtained for the orbital 
and spectroscopic parameters are in good agreement with those obtained by others using more 
traditional methods. 
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Figure 2: Radial velocities of the primary and secondary component of V373 Cas plotted against 
the photometric phase. The observational data are from Hill &: Fisher (1987). 
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Figure 3: Same as Fig. 2, for V2388 Oph. The observational data are from Rucinski et al. 
(2002). 
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Figure 4: Same as Fig. 2, for V401 Cyg. The observational data are from Rucinski et al. (2002). 
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Figure 5: Same as Fig. 2, for GM Dra. The observational data are from Rucinski et al. (2002). 



V523 Cas 




Photometric Phase 

Figure 6: Same as Fig. 2, for V523 Cas. The observational data are from Rucinski et al. (2003). 
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Figure 7: Same as Fig. 2, for AB And. The observational data are from Pych et al. (2004). 



HD 141929 
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Figure 8: Same as Fig. 2, for HD 141929. The observational data are from Carrier (2002). 
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Table 1: Orbital parameters of V373 Cas. 




This Paper 


Karami & Mohebi (2009) 


Hill & Fisher (1987) 


Primary 








V cm (kms~ v 


I -25 ±1 


-25.14 ±0.76 


-24.5 ± 2 


K p (kms^ 1 ) 


110 ± 1 


109.52 ±0.22 


106.7 ± 2.7 


e 


0.085 ±0.001 


0.0972 ± 0.0002 




o 


355 ±5 


344.5 ± 0.3 




Secondary 








V cm (kms^ v 


) -25 ±1 


-25.14 ±0.76 




K s (kms^ 1 ) 


146 ± 1 


145.53 ± 0.02 






C s — Cp 


€ s — Cp 




In 


175 ±5 


w s = w p - 180° 




Table 2: Combined spectroscopic elements of V373 Cas. 


Parameter 


This Paper 


Karami & Mohebi (2009) 


Hill & Fisher (1987) 


m p sin 3 %/Mq 


13.1 ±0.3 


12.98 ±0.03 


12.6 ±0.2 


m s sin 3 i/Mq 


9.9 ±0.2 


9.77 ±0.04 


9.3 ±0.2 


(a p + a s ) sini 


/Rq 67.5 ± 0.5 


67.3 ±0.1 


66.1 ±0.9 


m p /m s 


1.32 ±0.02 


1.33 ±0.03 


1.35 ±0.04 


Table 3: Same as Table 1, for V2388 Oph. 




This Paper 


Karami & Mohebi (2009) 


Rucinski et al. (2002) 


Primary 








V cm (kms^ 1 ) 


-25 ±1 


-25.35 ±0.88 


-25.88 ±0.52 


Kp (kms~ l ) 


45 ± 1 


44.71 ± 0.01 


44.62 ± 0.48 


e 


0.008 ± 0.001 


0.006 ± 0.001 




o 


285 ±5 


276 ±2 




Secondary 








V cm (kms^ 1 ) 


-25 ±1 


-25.35 ±0.88 


-25.88 ±0.52 


K s (kms^ 1 ) 


242 ± 1 


241.99 ± 0.03 


240.22 ± 0.98 




6s — 


CS = dp 




In 


105 ±5 


w s = w p - 180° 





Table 4: Same as Table 2, for V2388 Oph. 



Parameter This Paper Karami & Mohebi (2009) Rucinski et al. (2002) 



mpSin 3 i/M Q 


1.66 ± 0.03 


1.653 ±0.001 




m s sin 3 i/Mq 


0.31 ±0.01 


0.306 ± 0.001 




(a p + a s ) sin i/R Q 


4.55 ±0.03 


4.545 ± 0.001 




m s jm v 


0.187 ±0.005 


0.185 ±0.001 


0.186 ±0.002 


{m,p + m s ) sin 3 i/Mq 


1.97 ±0.04 


1.959 ± 0.002 


1.926 ± 0.030 
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Table 5: Same as Table 1, for V401 Cyg. 





This Paper 


Karami & Mohebi (2009) 


Rucinski et al. (2002) 


Primary 








V cm (kms^ 1 ) 


27 ± 1 


27.24 ±0.15 


25.53 ±2.14 


K p (kms^ 1 ) 


69 ±1 


68.46 ± 0.02 


72.23 ± 2.43 




0.006 ± 0.001 


0.0046 ± 0.0001 




In 


185 ±5 


181 ±5 




Secondary 








V cm (kms^ 1 ) 


27 ± 1 


27.24 ±0.15 


25.53 ±2.14 


K s (kms^ 1 ) 


248 ± 1 


247.99 ± 0.02 


249.13 ±4.53 




& s — &p 


e s — dp 




In 


5± 5 


w s = w p - 180° 





Table 6: Same as Table 2, for V401 Cyg. 



Parameter This Paper Karami & Mohebi (2009) Rucinski et al. (2002) 



m p sin 3 %/Mq 


1.50 ±0.02 


1.499 ± 0.001 




m s sin 3 i/Mq 


0.42 ± 0.01 


0.414 ± 0.001 




(dp + a s ) sin i/R® 


3.65 ± 0.02 


3.643 ± 0.001 




m s /m p 


0.278 ± 0.005 


0.276 ± 0.001 


0.290 ±0.011 


(m p + m s ) sin 3 i/M Q 


1.92 ±0.03 


1.913 ±0.002 


2.008 ±0.130 


Table 7: Same as Table 1, for GM Dra. 


This Paper Karami & Mohebi (2007a) Rucinski et al. (2002) 


Primary 








V cm (kms^ 1 ) 


10 ±1 


10.55 ±0.59 


9.12 ± 1.63 


Kp (kms~ l ) 


258 ± 1 


258.67 ±0.59 


258.72 ± 2.66 


e 0.005 ± 0.001 


0.008 ± 0.004 




O 


165 ±5 


175 ±2 




Secondary 








V cm (kms- 1 ) 


10 ±1 


10.55 ±0.59 


9.12 ± 1.63 


K s (kms^ 1 ) 


46 ±1 


46.59 ± 0.98 


46.69 ± 1.74 




€ s — Cp 


0.003 ±0.001 




o 


330 ±5 


339 ±2 




Table 8: Same as Table 2, for GM Dra. 


Parameter 


This Paper 


Karami & Mohebi (2007a) 


Rucinski et al. (2002) 


m p sin 3 i/Mq 


0.15 ±0.01 


0.15 ±0.01 




m s sin 3 i/Mq 


0.84 ±0.01 


0.85 ±0.01 




(dp + d s ) sin i/Rq 


2.03 ±0.01 


2.04 ±0.01 




m p / m s 


0.178 ±0.005 


0.176 ±0.004 


0.180 ±0.002 


{nip + ni s ) sin 3 i/Mq 


0.99 ± 0.02 


1.00 ±0.02 


1.002 ± 0.042 
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Same as Table 1, for V523 Cas. 




Ihis Paper 


Karami 6z Monebi (zUUyj 


Kucmski et al. (zOUoj 


Primary 








V cm {kms^ 1 ) 


-2±1 


-2.31 ±0.71 


-2.54 ±0.90 


K p (kms~ l ) 


237 ± 1 


236.22 ± 0.04 


235.95 ± 1.41 


e 


0.002 ± 0.001 


0.0012 ± 0.0002 




o 


190 ±5 


190 ± 11 




Secondary 








V cm (kms~ l ) 


-2±1 


-2.31 ±0.71 


-2.54 ±0.90 


K s (kms^ 1 ) 


123 ± 1 


122.38 ± 0.02 


121.64 ±1.14 


e 


e s — 6p 


— 6p 




o 


10 ±5 


w s = w p - 180° 





Table 10: Same as Table 2, for V523 Cas. 


Parameter 


This Paper 


Karami & Mohebi (2009) 


Rucinski et al. (2003) 


nip sin 3 i/Mq 


0.39 ±0.01 


0.381 ± 0.002 




m s sin 3 i/M@ 


0.74 ±0.01 


0.736 ±0.001 




(a p + a s ) smi/R Q 


1.66 ±0.01 


1.6557 ±0.0003 




m p /m s 


0.52 ±0.01 


0.5177 ±0.0002 


0.516 ±0.007 


{nip + m s ) sin 3 i/M Q 


1.13 ±0.02 


1.117 ±0.003 


1.11 ±0.24 



Table 11: Same as Table 1, for AB And. 



This Paper Karami et al. (2008) Pych et al. (2004) 



Primary 








V cm (kms~ l ) 


-27 ± 1 


-27.26 ±0.66 


-27.53 ±0.67 


Kp (kms~ l ) 


233 ±1 


232.69 ± 0.02 


232.88 ±0.83 




0.002 ± 0.001 


0.00109 ± 0.00005 




o 


220 ± 5 


231 ±3 




Secondary 








V cm (kms' 1 ) 


-27 ± 1 


-27.26 ±0.66 


-27.53 ±0.67 


K s (kms^ 1 ) 


133 ±1 


132.43 ±0.01 


130.32 ± 1.17 




6s = Gp 


C s — Gp 




In 


40 ±5 


LU s =LUp- 180° 





Table 12: Same as Table 2, for AB And. 



Parameter 


This Paper 


Karami et al. (2008) 


Pych et al. (2004) 


nip sin 3 i/M & 
m s sin 3 i/M & 
(a p + a s ) sin i/R Q 
nip/m s 

{nip + m s ) sin 3 i/M@ 


0.61 ±0.01 
1.07 ±0.02 
2.40 ± 0.01 
0.57 ±0.01 
1.68 ±0.03 


0.6071 ± 0.0001 
1.0668 ± 0.0002 
2.3943 ± 0.0002 
0.5691 ± 0.0001 
1.6739 ± 0.0003 


0.560 ± 0.007 
1.648 ±0.020 
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Table 13: Same as Table 1, for HD 141929. 




This Paper 


Karami & Mohebi (2007b) 


Carrier (2002) 


Primary 








V cm (kms- 1 ) 


-1 ± 1 


-0.44 ±0.12 


-0.33 ± 0.08 


K p (kms~ l ) 


11 ± 1 


10.38 ±0.01 


10.58 ±0.16 




0.398 ±0.001 


0.391 ± 0.001 


0.393 ± 0.008 


o 


140 ±5 


148.04 ± 0.34 


145.7 ± 1.7 


Secondary 








V cm (kms' 1 ) 


-1 ± 1 


-0.44 ±0.12 


-0.33 ± 0.08 


K s (kms~ l ) 


10 ±1 


9.87 ±0.01 


9.95 ±0.17 




€ s — €p 


0.389 ± 0.001 


0.393 ± 0.008 


In 


320 ±5 


324.03 ± 0.34 


325.7 ± 1.7 



Table 14: Same as Table 2, for HD 141929. 



Parameter This paper Karami & Mohebi (2007b) Carrier (2002) 



m p sin 3 i/Mq 


0.017 ±0.005 


0.0163 ± 0.0001 


0.01681 ± 0.00064 


m s sin 3 i/M & 


0.019 ±0.005 


0.0171 ±0.0001 


0.01789 ± 0.00067 




6.9 ±0.6 


6.53 ±0.01 


6.65 ±0.11 


a s sini/10 6 Km 


6.3 ±0.6 


6.21 ±0.02 


6.25 ±0.11 


m p /m s 


0.9 ±0.2 


0.953 ± 0.004 


0.94 ± 0.02 
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